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ABSTRACT 



' Context. We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact 

H II region (UCH ii) Mon R2. 

Aims. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate 
the chemistry of these species. 

Methods. We present an observational study of the small hydrocarbons CH, CCH and C-C3H2 in Mon R2 combining spectral mapping 
^ , data obtained with the IRAM-30m telescope and the Herschel space observatory. We determine the column densities of these species, 

and compare their spatial distributions with that of poly cyclic aromatic hydrocarbon (PAH), which trace the PDR. We compare the 
observational results with different chemical models to explore the relative importance of gas-phase, grain- surface and time-dependent 
chemistry in these environments. 

^ , Results. The emission of the small hydrocarbons show different spatial patterns. The CCH emission is extended while CH and C-C3H2 

-^j- . are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of C-C3H2 and CCH shows spatial 

00 ' variations up to a factor of a few, increasing from A/^(c-C3H2)/A/^(CCH) ^ 0.004 in the envelope to a maximum of ^ 0.015 - 0.029 

' towards the 8 yt/m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite 

, constant over a wide range of Go, whereas the abundance of C-C3H2 is higher in low-UV PDRs, with the A/^(c-C3H2)/A/^(CCH) ratio 

ranging ^^0.008-0.08 from high to low UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the 
abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce C-C3H2. In the low- 
density molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons 
! abundances. 

Conclusions. Our study shows that the small hydrocarbons CCH and C-C3H2 present a complex chemistry in which UV photons, 
grain-surface chemistry and time dependent effects contribute to determine their abundances. Each of these effects may be dominant 
depending on the local physical conditions, and the superposition of different regions along the line of sight leads to the variety of 
! measured abundances. 

^ . 

^ , Key words. ISM: abundances - ISM: individual objects: Mon R2 - Photon-dominated regions - ISM: molecules - Radio lines: ISM 



1. Introduction 

Ultracompact H 11 regions (UCH 11) represent one of the earliest 
phases in the formation of a massive star. In these environments, 
the chemistry is strongly driven by the mutual interaction be- 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 

Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



tween the gas, the dust and the UV photons from nearby stars. 
The H-ionizing photons are absorbed in a thin layer around the 
star, forming the H 11 region, while the UV radiation carrying en- 
ergies less than 13.6 eV penetrate in the deeper layers of the 
molecular cloud and produce the so-called photon-dominated 
region (PDR). In these environments, the thermal balance, the 
chemistry and ionization balance are driven by UV photons. 
PDRs associated with UCH 11 are characterized by extreme UV 
field intensiti es (expressed in units of the Habing field Go, see 
lHabing|[T968l) and small physical scales (< 0.1 pc). These envi- 
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ronments are therefore ideal to study the effect of a strong UV 
field on the chemistry of gas-phase species. 

The chemistry of small hydrocarbons in PDRs is still poorly 
constrained. The comparison of observations and gas-phase 
chemical models showed that this type of chemistry alone cannot 
account for the abundance of small hydrocarbons, especially for 
those species containing more than two C atoms: in NGC 7023, 
the Ori on Bar (Fuente et al. 2003), IC63 and the Horsehead 
nebula (iTeyssier et al.l |2004|) the abundances of species such 
as C-C3H2 and C4H derived from the observations are at least 
an order of magnitude larger compared to model predictions. 
This discrepancy was confirmed by high spatial resolutio n ob- 
serva tions and detailed modeling of the Horsehead nebula (IPetyl 
I2005h . in which the emission of small hydrocarbons was com- 
pared to other FDR tracers such as the mid-IR emission in the 
aromatic infrared bands (AIBs) due to poly cyclic aromatic hy- 
drocarbons (PAHs). The authors suggested that small hydrocar- 
bons may form through non gas-phase processes such a s the 
photo -destruction of the AIB carriers. In the Horsehead (iPetyl 
I2OO5I) and in NGC 7023 (Pilleri et al, in prep.), the spatial dis- 
tributions of CCH and C-C3H2 are strikingly similar. Similarly, 
a tight cor relation betwee n these two hydroc arbons has been ob- 
served bv lLucas&Lisztl (I2OO O) and Ge rin et al.l (l201lb on sev- 
eral lines of sights in the diff'use medium. The small radical 
CH has also a relatively constant abundance in these environ- 
ments ( Sheff'er et al. 2008). In low- to mild-UV irradiated PDRs 
(Go ~ 100 for the Horsehead and ~ 2600 for NGC7023) the 
tight correlation between CCH and C-C3H2 suggests a common 
origin of these two small hydrocarbons. However, this spatial 
correlation has not been explored yet in highly UV-illuminated 
PDRs (Go > 1 X 10^). 

In this paper, we present new observations of several lines of 
CH, CCH and C-C3H2 in the PDRs associated with the UCHii 
region Mon R2 and investigate observationally the relative vari- 
ations of these species. In Sect.[2lwe present the source and the 
observations, and the results are reported in Sect.|3] In Sect.lH 
we calculate column densities and abundances, while in Sect.[5] 
we compare the spatial distributions to that of other PDR tracers 
and discuss the quantitative consistence with diff'erent chemical 
models. Finally, Sect.[6] presents the conclusions and perspec- 
tives. 



2. Observations 

2.1. MonR2 

Mon R2 is a close-by ultracompact H 11 region (d=830 pc), which 
comprises several PDRs that can be spatially resolved in both the 
mm and IR domains. Due to its brightness and proximity, this 
source can be used as reference for other PDRs illuminated by a 
strong UV field. The most intense UV source is called IRS 1 and 
it is located at the center of the UCH 11 (see Fig[T]). A low-density 
molecular cloud surrounds this innermost region and extends for 
several arcminutes. The brightest PDR of Mon R2, which sur- 
rounds the UCH II and peaks at about 20'' to the NW of IRSl, 
shows both very intense AIB emission and H2 rotational lines 
dBerne et al. 2009) and is due to the illumination of the internal 
walls of the molecular cloud by IRS 1 . This PDR is very com- 
pact (d ~ 5'') and illuminated by a very strong radiation field 
(Go ~ 5 X 10^). A secondary PDR, detected north of the UCH 11 
at a distance of about 1 ' from IRS 1 , is much fainter and more ex- 
tended. The northern Molecular Peak posi tion at [0'' ;40''] (he re- 
after MP2 following the nomenclature of lGinard et al.ll2012l) IS 
clearly detected in the AIB emission at 8 /im, and its chemical 
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Fig. 1. Overview of th e region. The color image displays the 
IRAC 8yum emission (iGinard et al.l 12012 ^. in which both the 
bright circular PDR and the extended northern PDR can be dis- 
tinguished beneath the peaks due to the young stellar objects. 
The white contours represent the integrated intensity of the C^^O 
(2^ 1) lin e between 5 and 15 k m s"^ observed at the IRAM 30m 
telescope (iPilleri et al.l 12012a'), with contour levels from 16 to 
24 K km s"^ in linear steps of 2 K km s"^ The black squares rep- 
resent the positions of the ionization front (IF) and the molecular 
peak (MP2), whereas the triangles and the white line indicate the 
positions observed with HIFI. The black stars indicate the posi- 
tion of th e most intens e infrared sources following the nomen- 
clature of iHenning et al. (1992) . The ye llow arrow indicate s the 
direction of the NW molecular outflow (iTafalla etal.lll994 . 



properties see m to be similar to those of low- to mild-UV irra- 
diated PDRs ( Ginardetal.ll2012h . The origin of this PDR could 
be the outer envelope of the molecular cloud illuminated by an 
external UV field, or it could be due to UV photons that are es- 
caping from IRSl through a hole, illuminating its walls. 

Although there are some asymmetries in the general shape 
of the cloud, such as the molecular hole to the NW, we can as- 
sume as a zero-order approximation that the dense core around 
the UCH II region is spherical and composed by concentric shells 
centered on IRS 1 , a simplified model which is close to reality in 
most outwards directions. Each of the shells has its own physi- 
cal conditions (kinetic temperature, gas density. Go) which vary 
with radius. In the followi ng, we will assume the physical pa- 
rameters presented in Pilleri et al.l (l2012a ) and schematized in 
Fig.O the region immediately surrounding the star is an H 11 re- 
gion, which is expanding at a velocity of ~ 10 km s"^ (IChoi et al.l 
2000) and which is free of any molecular gas. The diameter of 
the region is ~ 40'', corresponding to ~ 0.08 pc. This is fol- 
lowed by a thin (1 mag), high UV-irradiated PDR, which has a 
density o f nu^ ~ 2x10^ cm ~^ and a gas kinetic temperature of 
~ 100 K jBeme et al.l l2009). Going farther from IRSl, there is 
a thicker (~ 8 m ag), warm (r ^in ^ 70 K) and high density (nu2 
~ 3 X 10^ cm-^ iRizzo et al.ll2QQ3l) PDR. which is expanding a t 
a velocity of < 1 km s'^ (iFuent e et al.ll2010l:lPilleri et al.ll2012ah . 
This high-density layer is surrounded by the parent molecular 
cloud, which is about 40 mag thick, has a relatively low density 
(nu2 ~ 5 X 10"^ cm~^) and temperature (T^kin ^ 35 K ), and is rather 
quiescent dFuente et al]l2010l: (pmeri et al.ll2012i) . The cloud is 
also illuminated from the outside with a UV field of Go ~ 100. 
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The presence of these different phases and its proximity 
make Mon R2 the best source to investigate the chemistry 
and physics of these extreme PDRs and their parent molecular 
clouds. 



observed in single pointing observing mode towards the IF, with 
a double beam switching pattern. In this paper we also compare 
the new obs ervation s with the [C 11] strip already presented by 
|Pillerietal.l(l2Q12ah . 



radial distance from IRS I 



HII region 

Molecular 'hole' 
r~0.08 pc 



Exposed PDR 

Go- 4 10^ 
I n(H2)~-2 10=cm-3 

Ar~ 0.001 pc 
(N(H2)~ 1x10^1 cm-- 
Tkm> 100 K 



High density PDR 

n(H2)~3 10^ cm- 

Ar~ 0.001 pc 
NCH2)= 1 10^^ cm-^ 
Tkm ~ 70 K 



Envelope 

n(H2)~5 lO^cm-^ 
Ar~0.26 pc 
N(H2)=5 10^2 
Tkm ~ 35 K 



Fig. 2. A ID schematic of the geometry of MonR2 as a function 
of the radial distance to IRS 1 . 



2.2. IRAM and Herschel observations 

We observed several low-lying rotational transitions of CCH and 
C-C3H2 at mm wavelengths using the IRAM-3Qm telesco pe. We 
completed this dataset with Herschel (iPilbratt et al .1120101) obser- 
vations of a higher excitation line of CCH and two lines of CH. 
A summary of the spectroscopic parameters for all the observed 
lines are given in Table IB. II Throughout the paper, we will use 
main beam temperature as intensity scale. Offset positions are 
given relative to the ionization front (IF: Qf720oo=06h07m46.2s, 
<^/20oo = -06°23'08.3'0. For all the observations, the OFF posi- 
tion was chosen to be free of molecular emission ([-i-10',0'] for 
Herschel; [+400'', -400''] for IRAM). A summary of the obser- 
vations is shown in Table [T] 

The 3mm and 2mm observations were performed in July 
2012 at the IRAM-30m telescope in Pico Veleta (Spain). We 
observed the CCH multiplet at 87 GHz and the C-C3H2 line 
at 85 GHz using the EMIR receivers with the VESPA corre- 
lators, covering a 180"xl80" region centered toward the IF. 
The 1 mm observations were performed in two different ses- 
sions using the HERA 3x3 dual-polarisation receiver and in a 
successive run with the EMIR receivers. During the first run 
(HERA) we obtained a 150"xl50" map of the CCH multiplet 
at 262 GHz and the C-C3H2 line at 217 GHz. The second ob- 
serving run (January 201 2) was part o f a 2-D spectral survey 
at 1mm (Trevino-Morales, et al.ll2012l in prep.). The resulting 
C-C3H2 maps have a lower signal-to-noise ratio (SNR) and the 
field of view covered in this second run was smaller, 2'x2'. 

The Herschel observations were obtained with the 
Heterodyne Instrum ent for the Far Infrared (HIFI, 
Ide Graauw et al.l 12010) in the context of the W ADI guaranteed- 
time key program (lOssenkopf et al.l l201lb . The CH lines 
were observed with both the wide band spectrometer (WBS) 
and the high-resolution spectrometer (HRS). They provide 
spectral resolutions (at 500 GHz) of Av = 1.2 km s"^ and 
0.3kms"\ respectively, allowing to resolve the line profiles. 
Cross-calibration between the two backends and between the 
two polarizations indicates an uncertainty of ~ 20%. The CH 
(J = 3/2-1/2) and CCH (N=6-5) observations were performed 
in on-the-fly (OTF) observing mode along the strip indicated 
in Fig. [T] while the excited (/ = 5/2 - 3/2) line of CH was 



2.3. Data reduction 

The basic data reduction of Herschel observations was per- 
formed us ing the st andard pipeline provided with the version 7.0 
of HIP^l (IOttll2010 ) and Level 2 data were exported to the FITS 
format. The Herschel and IRAM-30m data were then processed 
using the GILDAS/CLASS softwar^ suite (iPetvll200 5). The spec- 
tra were first calibrated to the T* scale using the chopper wheel 
method (iPenzias & Burruslll973l) . and finally converted to the 
main beam temperatures (r^b) using the nominal forward (Feff ) 
and main beam (^eff) efficiencies (in Tabled we show the value 
r]i = Bq^/Fq^ ). The HIFI efficie ncies have been recently calcu- 
lated by .Roelfsema et al.l (l2012l) . The resulting amplitude accu- 
racy is ~ 10% for both HIFI and IRAM observations. 

Visual inspection of the IRAM data revealed for all the ob- 
served lines 1) the presence of more or less pronounced spikes 
at one quarter and three quarter of the correlator window, 2) 
the presence of platforming in the middle of the correlator win- 
dow, 3) but otherwise clean baselines. Windows around the well- 
defined frequencies displaying spikes were set before baselin- 
ing. The line windows were defined on the spectra averaged 
over the map, from [-1,24.5 km s"^] for the brightest lines to 
[6,16 km s"^] for the faintest ones. To correct for the platform- 
ing, a zero-order baseline was subtracted from each correlator 
sub-bands of each spectrum. A first-order baseline was needed 
for C-C3H2. The bad channels, which occurred in a line-free sec- 
tion of the HERA spectra, were replaced by Gaussian noise of 
same rms. The resulting spectra were resampled and convolved 
with a Gaussian to a map resolution larger than the HPBW of the 
telescope. The map resolution is reported in Table [T] The spectral 
cubes were finally blanked where the spectra weights (resulting 
from the gridding) were too low, giving the fringed aspect of the 
map edges. 



3. Results 

In this section, we describe the spectral and spatial variations 
of the various molecular tracers. First, we present all the obser- 
vations obtained towards the IF, whic h is one of the best-known 
position (see, for in stance, iFuente et a l. 2010; Ginard et al .120121: 
Pilleri et al.l l2012al) . We then describe the spatial/spectral varia- 
tions along the Herschel strip (see Fig.O and finally, we present 
the full maps. 



3.1. The ionization front 

In Fig. [3] we show the spectra of several CH, CCH and C-C3H2 
transitions observed towards the IF with the IRAM-30m and 
Herschel. The lower-lying transitions show a self- absorption 
feature at the ambient velocity of the cloud (~ 11 kms"^). It is 
clearly detected in the main hyperfine components of CH, CCH 
and C-C3H2, and coincides with the se lf- absorption detec ted in 
the [Cii] strip observed with Herschel (iPilleri et al.ll2012al) . The 



^ RIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia 

^ See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILD AS softwares. 
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Table 1. Summary of the observations 



Line(s) 






Frequency 


Instrument 


FWHM 




Position 


OM^ 


T c 

rms 












L J 










C^^O / = 2 - 


1 




219.560 


IRAM-HERA 


10.0 


0.67 


Map 


OTF 


0.025 


C^^O 7 = 5- 


4 




548.830 


Herschel-WFl 


38.6 


0.79 


Strip 


OTF 


0.008 


CH7 = 3/2- 


1/2 




536.761 * 


Herschel-mFl 


38 


0.77 


Strip 


OTF 


0.041 


CH 7 = 5/2 - 


3/2 




1656.961 * 


Herschel-mVl 


12.8 


0.74 


IF 


SP 


0.284 


CCH = 1 - 







87.317 * 


IRAM 


29 


0.85 


Map 


OTF 


0.070 


CCH = 3 - 


2 




262.004 * 


IRAM-HERA 


9 


0.6 


Map 


OTF 


0.380 


v^v^n. iV — u 


5 




523.972 * 




40.5 


0.78 


Strip 


OTF 


0.057 


C-C3H2 Jka,Kc 


= 2i,2 


-lo,i 


85.33889 


IRAM 


29 


0.85 


Map 


OTF 


0.070 


C-C3H2 Jka,Kc 


= 4o,4 


-3i,3 


150.82067 


IRAM 


16 


0.64 


Map 


OTF 


0.054 


C-C3H2 Jka,Kc 


= 4i,4 


- 3o,3 


150.85191 


IRAM 


16 


0.64 


Map 


OTF 


0.054 


C-C3H2 Jka,Kc 


= 61,6 


- 5o,5 


217.82215* 


IRAM-HERA 


11 


0.61 


Map 


OTF 


0.016 


C-C3H2 Jka,Kc 


- 5i,4 


-42,3 


217.94005 


IRAM 


11 


0.61 


Map 


OTF 


0.150 


C-C3H2 Jka,Kc 


= 43,2 


-42,1 


227.16913 


IRAM 


10 


0.59 


Map 


OTF 


0.150 


C-C3H2 Jka,Kc 




- 61,6 


251.31434* 


IRAM 


9 


0.55 


Map 


OTF 


0.115 


C-C3H2 Jka,Kc 


= 62,5 


- 5i,6 


251.52730 


IRAM 


9 


0.55 


Map 


OTF 


0.115 



T], = B,ff/F,ff. 

Observing Modes: on the fly (OTF), single pointing (SP). 
'^''^ Calculated on a scale with Av = 0.5 km s~^ 
*^*^ Frequency of the most intense transition of the hyperfine structure. 



self-absorption is barely visible in the faintest hyperfine com- 
ponents of the CCH 1-0 line and is not detected in any of the 
higher energy transitions. The structure of the CCH 3-2 line is 
due to the superposition of two hyperfine transitions with very 
similar frequencies that are not self-absorbedEl. The lower-lying 
transitions show some emission in the r ed and blue wings , which 
may trace either an expanding PD RjiPilleri et al.ir2Q12ah or the 
molecular outflow (.Tafalla et al... 19971) . 

3.2. The NE-SW strip 

CH and C^^O (5-4) have been observed only towards the NE-SW 
strip sh own in Fig.[T] This strip, already studied in iPilleri et aP 
(I2012ah . passes through the IF and extend up to 1.5' in each di- 
rection, crossing the main PDRs of the region and extending to 
ofl'sets where the emission is dominated by the molecular en- 
velope. However, the large beam of Herschel at this frequency 
(~ 4C) results in even the positions [+40'', ±40"] picking-up 
some emission from the innermost PDR. 

In Fig. 131 we show the position- velocity (PV) diagrams along 
the cut observed with Hersch el, i.e. in the direction per pendic- 
ular to the molecular outflow (iTafalla etal.lll994l Il997h . which 
extends in the northwest- southeast direction from the IF. For a 
comparison, we als o display the PV d iagram of the [C 11] line 
along the same cut (IPilleri et al.ll2012ah . For each transition, we 
show the diagram of the component with the highest signal-to- 
noise ratio (SNR). Although the data have very different angular 
resolutions, the spatial variation of CH is in some way similar 
to that of [C 11] : for instance, in both tracers the emission is ex- 
tended up to 60" from the IF and the main component shows 
self-absorption through the whole cut. Both lines present a sec- 



^ This is suggested by the excellent quality of the fit of the observed 
lines using a Gaussian function for each of the hyperfine component, 
and assuming a standard intensity ratio and frequency shift, see Table 



ondary peak at a distance of ~ 10" at about Skms"^. The line 
widths are also similar indicating that the CH emission, like 
[Cii], originates in the innermost layers of the expanding cen- 
tral PDR. 

The PV diagram of CCH 1-0 is very similar to the C^^O 
2-1 line in terms of the spatial extension and line width. They 
are both detected along the entire strip, and their typical line 
width is about 2.5 km s"^. However, the CCH line toward the IF 
displays self- absorption at central velocities and emission in the 
red wing. Similarly to the CCH 1-0 line, the C-C3H2 2 1,2 - lo,i 
transition shows self- absorption at 1 1 km s"^ and emission in the 
red wing toward the IF. However, it is less spatially extended 
and therefore likely associated with the innermost layers of the 
PDR. The higher-energy transitions of all species are less spa- 
tially extended compared to [C 11] and CH and do not show self- 
absorption. Finally, there is a systematic shift of the peak ve- 
locity of about 2kms"^ going from the south-west towards the 
north-east of the IF, which has been attributed to either an asym- 
metry in the expan sion pattern or to a large-scale s low rotation 
of the entire cloud (iLorenll 19771 IPilleri et al.ll2Q12ab . 

3.3. Maps 

The maps carried out with the 30m telescope allow us to inves- 
tigate about the spatial distribution of CCH and C-C3H2. Figure 
[5] shows the spatial distribution of the integrated line intensity of 
the CCH 3-2 and C-C3H2 6-5 lines between 5 and 15 km s"^ ob- 
served with HERA0 at the spatial resolution of 10". The spatial 
pattern of CCH is very similar to that of C^^O 2-1 (Fig.Q): the 
emission to the northeast and southwest extends up to V from 
the IF, and peaks at a distance of ~ 30". The emission is opened 
to the north and the northwest, following the cometary shape of 
the molecular bar. The C-C3H2 emission is more compact, with 
the center at an off'set [-10", +10"] relative to the IF. The emis- 



^ The EMIR maps of all the observed lines are shown in Fig. IB. II 
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Fig. 3. Herschel and IRAM observations of CCH and C-C3H2 towards the IF. 1mm observations have been convolved to the spatial 
resolution of 39''. Solid blue lines show the single-component Gaussian fits. The velocity scale in each panel is relative to the most 
intense line in the hyperfine groups (see also Tables [Tl and IB. II) . 



sion extends less than 40'' from the center, and peaks at a dis- 
tance of about 10" south from the IF. The line intensity also 
displays a peak toward the MP2 position, which does not have 
a correspondence in the CCH map. The emission to the north 

is more extende d compared to CCH, following the IRAC 8yum ^-^^ The SW-NE strip 



spatial distribution of these tracers. The self- absorption features 
are discussed in Appendix lAl 



emission dGinard et al.ll2012l) . 



4. Column densities and abundances 

In this section, we derive the column densities and abundances 
of CH, CCH and C-C3H2. First, we concentrate in the SW-NE 
strip, where we use the Herschel observations to study CH and 
the excitation of C^^O. We study in more details the excitation of 
CCH using a two-phase LVG model towards the IF and study the 
emission in the line wings. Finally, we extrapolate these results 
to the analysis of the full maps to obtain more information on the 



To obtain an estimate of the excitation temperature and column 
density of C^^O along the Herschel strip, we have combined 
our J=2-l o bservation with the J = 5-4 observations published 
in IPilleri et al.l (l2012ai) and with the J = 1-0 observations of 
iGinard et"aD (12012 ). We have convolved the data to the lowest 
spatial resolution (40" for the J = 5-4) and built rotational dia- 
grams for the five positions along the strip. The J=l-0 data can- 
not be convolved, since they were single pointing observations, 
but the spatial resolutions are not very diff'erent (29" vs 40") and 
the C^^O emission is very extended, so that the emission fills the 
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Fig. 4. PV diagrams of the Herschel and IRAM observations along the SW-NE strip. The Cii diagrams comes from iPilleri et al.l 
(l2Q12ah . The double peaked structure of the CCH lines is due to two hyperfine component of the same line, while the strong dip in 
both CH and [C 11] are partially due to self- absorption. The white line in the top-right corner of each diagram represents the HPBW 
of the observation. The black dotted line indicates vlsr = 11-5 km s"^. 



beam at both spatial resolutions. The corresponding rotational 
diagram s are shown i n Fig JB. 2 1 The results are consistent with 
those of iGinard et al.l (|2012): the rotation temperature slightly 
decreases from ~ 37 K toward the IF to ~ 31 K at a distance 
of 50''. We obtain a column density of 8.3 x 10^^ cm"^ toward 
the IF, and slightly (20%) higher values at larger offsets. These 
values are consistent within a factor of 2 with the geometrical 
model presented by Pilleri et al. ( 2012a) to fit the high- J lines of 
CO, ^-^CO and C^^O. This model consists of a spherical cloud 
with an Ay = 8 PDR and an Ay = 40 low-density envelope. This 



geometrical model yields a factor of 2 higher column densities 
compared to the C^^O values, assuming an isotopic ratio ^^O/^^O 
= 500 and a standard abundance of CO relative to H2 of 1 x 10""^. 
In the following, we assume therefore an uncertainty of a factor 
of 2 in the column density of H2 and in the abundances of the 
hydrocarbons. However, the absolute value of the column densi- 
ties, and consequently the relative abundances, are not aff'ected 
by this uncertainty. 
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studied in this paper. White contours in the left panel represent the C-C3H2 line integrated intensities (in the range 5-15 km s ^) 
starting from 0.6 km s"^ and in linear steps of 0.2 km s"^. 



Table 2. Integrated intensity (using scale) between 5 and 15 kms"\ after smoothing the data to the larger HPBW (W for 
CCH and C^^O, 28'' for C-C3H2). Non-detections are given at a 3cr level assuming a line width of 1 km s"^. See Table IB. II for the 
spectroscopic parameters. 





IF [0";0"] 


[20"; 20"] 


[40"; 40"] 


[-20"; -20"] 


[_40"; -40"] 


MP2 [0";40"] 


Transition 


Area 


Area 


Area 


Area 


Area 


Area 




[K kms"^] 


[K kms-^] 


[K kms-i] 


[K kms-^] 


[K kms"^] 


[K kms-^] 


Ci^0 2-1 


14.56 


15.55 


15.71 


15.68 


13.12 


13.90 


C^^O 5-4 


13.04 


13.38 


13.26 


12.12 


11.03 




CH 7=3/2- 1/2,/'=!- -0+ [5-12]*^^^ 


1.81 


1.43 


0.62 


1.03 


0.22 




CH, =3/2-1/2,^=2-1- [12-15]^^) 


1.19 


0.74 


0.36 


0.95 


0.22 




CH 7=5/2-3/2 


<0.84 












CCH ^=i_o,i7=o-i 


0.98 


1.35 


1.39 


1.61 


1.87 


1.29 


CCH N=3-2,Ju=5/2,F=4-3&F=3-2 


13.75 


13.85 


11.70 


15.63 


11.99 


12.19 


CCH A^=6-5,y„ = 13/2,F=7-6&F=6-5 


4.97 


4.97 


4.05 


5.02 


3.73 




C-C3H2 212 


0.73 


1.55 


1.71 


1.73 


1.32 


1.71 


C-C3H2 404 


0.45 


0.44 


0.44 


0.45 


0.19 


0.67 


C-C3H2 414 


1.42 


1.71 


1.34 


1.52 


0.80 


2.0 


C-C3H2 432 


0.63 


0.46 


0.37 


0.49 


0.32 


0.53 


C-C3H2 514 


0.76 


0.56 


0.33 


0.46 


0.20 


0.87 


C-C3H2 616 & 606 


0.99 


1.00 


0.52 


0.85 


0.42 


1.01 


C-C3H2 625 


0.57 


<0.3 


< 0.3 


<0.3 


< 0.3 


0.41 


C-C3H2 707 & 717 


0.69 


0.61 


< 0.4 


0.58 


<0.3 


0.78 



(a) jYiQ numbers in brackets indicate the velocity range over which the line has been integrated to avoid self-absorption or contamination from a 
satellite line. 



CH 

As shown in Fig.[3l the main component of CH is self-absorbed 
at the ambient cloud velocity and the line wings of the different 
components overlap. To obtain the integrated intensity of the CH 
line, we used the P = 1" - 0^ and the P = 2" - 1^ transitions 
for the velocity intervals [5-12] km s"^ and [12-15] km s"\ re- 
spectively. Using the first satellite to obtain the area in the blue 
interval we avoid the self- absorption that aff'ects the main com- 
ponent. On the other hand, the red part of the P = 1" - 0^ line 



is contaminated by the overlap with the second satellite, so that 
for this interval we have to rely on the main component. The 
integrated areas for each interval are shown in Table |2] 

When computing the CH column density, we face the diffi- 
culty that there are no collisional rates available for CH, so that 
one can only assume a given rotational temperature T^ot per ve- 
locity interval. The non-detection of the CH triplet at 1656 GHz 
provides an upper limit to T^ot if we assume LTE conditions, so 
that both transitions share the same rotation temperature. This 
gives upper limits for T^ot of 25 K (blue) and 30 K (red) towards 
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the IF. Considering the high dipole moment of CH and the phys- 
ical conditions in Mon R2, we can estimate a lower limit to the 
rotational temperature to be about 7 K. Assuming that the emis- 
sion is optically thin and that it fills the beam, we derive a column 
density of 1.2+^-^x10^4 cj^-2^Q^^j,^g^j^^jpfQj,j^^^^ ^ 10K,with 

the uncertainties driven by upper and lower limit in Trot- At the 
[-40''; -40''] position, the CH column density is a factor of 10 
lower, 1.7^^1 x 10^^ cm"^, assuming Trot = 10 K. At this posi- 
tion, Trot is likely lower compared to the IF, due to a decrease in 
kinetic temperature and local density, which yields an overesti- 
mate of the column density by a factor of 2 for Trot = 7 K. At 
off'sets larger than 60", the ground state transition of CH is not 
detected, implying that the column density is even lower. Thus, 
we consider the column density at [-40"; -40"] as an upper limit 
of the CH abundance in the envelope. 

Using the H column density derived through C^^O (N(tl) = 
9x10^^ cm"^), we estimate a CH abundance relative to H nuclei 
of ~ 1.5 X 10"^ toward the IF decreasing to ~ 3 x 10"^^ at an 
off'set of [-40", -40"], for T^ot = 10 K. This is consistent with 
most of the CH emission at central off'sets arising from the cen- 
tral expanding FDR. However, if the bulk of CH emission arises 
indeed in the PDR layers, using C^^O as a probe of H2 most 
likely underestimates the abundance because /) the abundance of 
C^^O is aff'ected by photodissociation and //) because C^^O may 
not coexist with the CH present in the gas where C^^O is partly 
photo-dissociated. If we consider that the emission stems only 
from the innermost 10 mag in each side of the spherical clou d 
(A^(H) = 4 X 10^2 Sect. I4:2l[5:2l and Pilleri et al.ll2012al) . 

the mean abundance in the PDR would be of 3. 3 ^ 10"^ relsi- 
tive to H nuclei, with the error bars driven by the uncertainty in 

^rot- 

CCH 

The CCH line at 87 GHz is composed of 6 hyperfine compo- 
nents. The most intense line is self-absorbed at a velocity of 
-11 km s"^. To estimate the molecular column density, we used 
the faintest observed component (87.407 GHz), which is less af- 
fected by self-absorption. We have observed four hyperfine com- 
ponents of the J = 3-2 and 6-5 multiplets, which are blended in 
pairs of two. For each group, we derived the total integrated in- 
tensity of the lines and then scaled it according to the relative 
strengths to estimate the contribution of each blended lines sep- 
arately. 

Figure lR2l shows the rotational diagrams at several position 
along the Herschel cut. We have convolved the J = 3-2 line to 
the spatial resolution of the 6-5 observations (39"). For all the 
lines, we assumed a 30% uncertainty for the intensity, except for 
the 1-0 line where we have assumed an uncertainty of a factor 
of 2 due to a possible self-absorption. For the 3mm transitions 
we used the hyperfine fit task inside CLASS, and found that the 
relative line intensities of the 3mm transitions are consistent with 
optically thin emission. For the higher frequencies, the hyperfine 
lines overlap and this method is not reliable. However, our LVG 
results (see below) show that the opacities of these lines are also 
lower than 0.1. 

One single component cannot fit the rotational diagram at 
any of the positions. For this reason, we used a two-components 
rotational diagrams to fit the integrated intensities of all the lines 
along the strip. The warm component has a rotational tempera- 
ture of ~ 35 K while the cold component is at about 10 K. These 
temperatures are consistent with the superposition of two phases 
along the line of sight: a dense an d high-UV i lluminated PDR 
and a low-UV irradiated envelope ( Pilleri et al.ll2012a) . 



The CCH column density is approximately constant (~ 7.5 x 
10^"^ cm"^) along the entire strip, but the column density of the 
hot component decreases by about 25% at larger off'sets, where 
the emission is dominated by the envelope. This supports the in- 
terpretation of the warm component being associated with the 
PDR around the H 11 region (hereafter, we will refer to this com- 
ponent as the PDR component). It is important to note that this 
gradient is a lower limit to the real one because these column 
densities have been determined at a low spatial resolution (40"), 
thus diluting higher gradients in their distribution that may be 
present in the beam. Using the H2 column density derived from 
C^^O, we estimate a mean CCH abundance of 9 x 10"^ relative 
to H nuclei. 

C-C3H2 

The C-C3H2 2 1,2 - 1 0,1 transition is self-absorbed towards the IF 
implying that we can only establish a lower limit to the column 
density of C-C3H2 at this position. Elsewhere, the spectra are not 
self-absorbed and we can derive accurate column densities. 

We convolved the higher- J transition to the spatial resolution 
of the 2 1,2 - 1 0,1 line, and used the rotational diagram technique 
to compute the corresponding column densities. The rotational 
diagrams show that the column density decreases at larger off'- 
sets, from ~ 6.0 x 10^^ cm"^ towards the IF to ~ 4.3 x 10^^ cm"^ 
at an off'set of [-40",-40"] to the SW. Only one component is 
needed to fit reasonably well the observed intensities, with a ro- 
tational temperature in the range 8-16 K. The excitation temper- 
ature seems to decrease farther from the IF. The spatial distribu- 
tion of this species suggest that its emission is dominated by the 
dense PDR around the H 11 region, similarly to the warm compo- 
nent of CCH. The diff'erence between the rotation temperature 
of C-C3H2 and that of the hot component of CCH is partially due 
to the higher diff'erent dipole moment of C-C3H2, which makes it 
much more difficult to excite by collisions. This is also consis- 
tent with the results of the LVG modeling in Sect. 14.21 

Again, using C^^O as a tracer of the H2 column density 
yields an abundance relative to H of X(c-C3H2) = 7.3 x 
10"^^ towards the IF, decreasing by about 30% at the position 
[-40", -40"]. As for CCH, more details can be obtained assum- 
ing a rotational temperature and using only the 1mm observa- 
tions, improving the spatial resolution by a factor of 2.5. Such 
improvement will be shown in Sect. l4.4l 

4.2. Two-phase decomposition towards the IF 

The CCH rotational diagrams indicate that there are two emis- 
sion components towards t he IF. This i s also co nsistent with the 
source model proposed by lPilleri et al.l (l2012al) . in which a very 
dense PDR (nu2 = 3x10^ cm"^) is surrounded by a lower density 
envelope (T^in = 35 K, =5x10^ cm"^). In the previous sec- 
tion we interpreted the warm component of CCH and the emis- 
sion of C-C3H2 as arising in the PDR, with the cold component 
arising from the envelope. In order to verify this assumption and 
quantify the excitation eff'ects on CCH and C-C3H2, we used the 
MADEX large velocity gradient (LVG) code ( Cemicharo 20T2I) 
to reproduce the CCH and C-C3H2 emission toward s the I F us- 
ing the collisional coeff icients of Spielfiedel et all (1201 2h and 
IChandra & Kegell (l2000l) . respectively. 

We only have a small number of transitions with diff'erent up- 
per energy levels, and therefore we cannot accurately constrain 
all the physical parameters (nu2, T]^n, N). Therefore , we have to 
rely on previous knowledge of the source structure dPilleri et al.l 
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l2Q12ah to fix the temperature and density in each of these two 
phases and fit the column densities. We assumed densities of 
=3x10^ cm"^ and 5x10"^ cm"^ and kinetic temperatures of 70 K 
and 35 K for the PDR and the envelope, respectively. The line 
width for both components has been set to 2 km s"^. By assum- 
ing these parameters and varying the column densities in each 
of the two phases, we obtain the best fit to the observations with 
the solutions shown in Table O The line opacities returned by 
MADEX are lower than 0.1 for all the transitions. For CCH, we 
reproduce all the observed intensities within a 10% uncertainty. 
For C-C3H2, we are able to reproduce within a 15% error the 
observations with the best SNR, while the other transitions are 
fitted within 30%. This approach gives the column densities and 
abundances reported in Table [3] The results of MADEX mod- 
eling for all the transitions observed toward the IF are shown 
on the rotational diagrams in Fig. IB. 21 The uncertainties in the 
abundances are dominated again by the uncertainty in the total 
H column density, and can be estimated to be a factor of 2. 

Table 3. Results of the two-phase LVG modeling towards the IF. 
Abundances are relative to H nuclei and are calculated assuming 
that the PDR has a thickness of Ay = 8 and the envelope of 
Ay = 40 on each side of the spherical cloud. 





Combined 


PDR 


Envelope 


N(CCH) [cm-2] 


1.0 X lO^^ 


5.0 X 10^4 


5.0 X 1014 


X(CCH) 


5.7 X 10-9 


1.7 X 10-^ 


3.5 X 10-9 


A^(c-C3H2)[cm-2] 


6.0 X 10^2 


4.0 X 10^2 


2.0 X 1012 


X(C-C3H2) 


3.5 X 10-11 


1.4 X 10-10 


1.4 X 10-11 


A^(c-C3H2)MCCH) 


0.006 


0.008 


0.004 



4.3. Emission in the high-velocity wings 

The ground- state transitions of both CCH and C-C3H2 show a 
broad emission at red-shifted velocities, which is not clearly 
detected in the higher-energy lines (see Fig. [3]). Toward the 
IF, the integrated intensity in the interval [15-25] km s"^ is 
0.64Kkms-i and 0.44 Kkms"^ for CCH and C-C3H2, respec- 
tively. This allows to estimate an upper limit to the gas density 
of the region emitting in the wings from LVG models that fit the 
intensity of the 3 mm lines and stay below the detection limit for 
the other transitions. This yields an upper limit to the gas density 
of = 1 X 10^ cm"^ for Tkin = 70 K. Assuming this density, 
we derive column densities of 3.5 x 10^^ cm"^ and 1.5 x 10^^ 
cm"^ for CCH and C-C3H2, respectively. 

We repeated the rotational diagram for the C^^O lines in this 
velocity interval to estimate the column density of the gas in this 
region, finding N(C^^O) = 1.9 x 10^"^ cm"^, which corresponds 
to N(B) = 1.9 X 10^^ cm"^, or about 1 mag. This can be used to 
obtain an estimate of the abundance of CCH and C-C3H2 in the 
region emitting in the wings, which gives 1.8 x 10"^ for CCH 
and 7.8 x 10"^^ for C-C3H2. If the density was lower (e.g. = 
1x10"^ cm"^), the column density required to reproduce the wing 
intensity would be 15 % higher for CCH and a factor of 4 higher 
for C-C3H2, yielding a higher A^(c-C3H2)MCCH) ratio. 

4.4. Spatial distributions 

In the previous paragraphs we studied the excitation of single 
positions for all the species and derived some patterns in both the 



column density and rotational temperature. Here, we generalize 
the study to the mapping observations. First, we convolved all 
the maps to the lowest spatial resolution in the mapping dataset 
(29'0- For each pixel in the map, we constructed a rotational 
diagram for each molecule and derived an excitation temperature 
and a column density for both CCH and C-C3H2. For simplicity, 
we used only one slope also for CCH. The results are shown in 
Fig. [61 

For C^^O, our analysis is limited by the fact that the J = 5-4 
line was observed only in the NE-SW strip. We have therefore 
to assume a given T^ot in each pixel based on this dataset. We 
have interpolated the mean Trot obtained along the strip to obtain 
the rotation temperature of C^^O as a function of the distance 
from the IF. Then we derived the corresponding column densities 
(Fig.O assuming LTE. Since the C^^O emission is optically thin, 
this yields an estimate of the H2 column density and allows to 
normalize the column densities obtained for CCH and C-C3H2 to 
derive their abundances. The corresponding maps are shown in 
Fig.[6l 

To have a deeper insight into the hydrocarbon chemistry 
we produced a map of A/^(c-C3H2)/A/^(CCH), the ratio of the 
C-C3H2 to CCH column densities. The A^(c-C3H2)MCCH) ra- 
tio shows variation up to a factor of 3 across the whole region, 
being minimum in the bulk of the molecular cloud, then in- 
creasing towards the PDR around the HII region and with the 
maximum located 40'' towards the North, i.e., at t he position of 
the low-UV PDR identified bv lGinard et al.1 (12012) . In fact, the 
A/^(c-C3H2)/A/^(CCH) seems to follow the emission of the hot dust 
(PAH and larger grains) traced by the 8 yum emission. We can im- 
prove the angular resolution of our A/^(c-C3H2)/A/^(CCH) image 
by using directly the CCH 3-2 and C-C3H2 6-5 images (both with 
an angular resolution of -10") to derive the A^(c-C3H2)/A/^(CCH) 
but assuming the CCH and C-C3H2 rotation temperatures shown 
in Fig. [6] although these rotation temperatures were derived with 
a lower (-29") angular resolution. Note that the angular reso- 
lution of the IRAC 8yum image is T\ This higher resolution 
A/^(c-C3H2)/A/^(CCH) image confirms the correlation with the 
8yum image (see Fig. [7]) corroborating the association of an en- 
hanced A^(c-C 3H2)MCCH) with the external part of the PDR. 
In Sect. 15.2.11 we will discuss the chemical implications of this 
correlation. 



5. Discussion 

5.1. Comparison with other PDRs 

In the previous section, we have derived the abundances of small 
hydrocarbons in the diff'erent environments of Mon R2. In the 
PDR, the abundances are X(CH) = 3.3 x 10"^ X(CCH) = 1.7 x 
10-^ andX(c-C3H2) = 1.4 x IQ-^^. 

Table |4] shows the comparison of the CCH and C-C3H2 abun- 
dances in diff'erent PDRs. With the exception of the Orion Bar, 
the CCH abundance is relatively constant in PDRs, in the range 
(0.6 - 1.4) X 10"^ The abundance in the PDR of Mon R2 
is consistent with these results. However, the en velope better 
matches with the values found in the Orion Bar (iFuente et al.l 
1996; Ungerechts et al. 1997; Fuente et al. 2003). The case of 
the Orion Bar is interesting because, similarly to Mon R2, it is 
a high-UV PDR with an associated envelope, the background 
molecular cloud not directly aff'ected by the UV photons from 
the Trapezium cluster. The abundances for the Orion Bar shown 
in Table m were estimated using low angular resolution observa- 
tions (29'')- The beam of these observations comprises both the 
high-UV irradiated PDR and the warm envelope which hinders 
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Fig. 6. First two rows: Maps of the integrated intensity of the CCH (262.004 GHz) and C-C3H2 (217.822 GHz) lines, and of the 
column density and rotation temperature derived applying the rotational diagram technique to each pixel. Black contours indicate 
the CCH 3-2 integra ted intensity, and go from 5 K km s"^ to 20 K km s~^ in linear st eps of 2.5 K km s"^ . The position of the infrared 
sources IRS 1/2/3/4 (iHenning et al.lll992 ) and the MP2 position iGinard et al.ll2012h are shown as black triangles. All the maps have 
been convolved to a spatial resolution of 29''. Bottom row: CCH and C-C3H2 abundance maps, obtained using C^^O as a reference 
and assuming ^^O/^^O = 500 and a CO abundance of 5 x 10"^ relative to H nuclei. The A/^(c-C3H2)/A/^(CCH) ratio is also shown. 



from a reliable comparison with the values in other PDRs. Note, 
however, that the CCH abundance in the Orion Bar is similar to 
the mean CCH abundance in Mon R2 (Table [S]). 



On the other side, the C-C3H2 abundance shows a larger vari- 
ation amongst the various PDRs. Whereas the PDRs with a rel- 
atively low Go show an abundance of (0.6 - 1.1) x 10"^, NGC 
7023 NW and Mon R2 show much lower values, (0.4 - 2) x 10"^^. 
This is also seen as the trend in the A/^(c-C3H2)/A/^(CCH) ratio in 
the different PDRs. Again, the values in the Orion Bar are more 
similar to the mean value in Mon R2. The A^(c-C3H2)/A^(CCH) 
ratio obtained in the PDRs of Mon R2 is an order of magnitude 
lower compared to that found in "classical" PDRs (s ee, for in- 
stance, Pety 2005) and in the diffuse medium (iGerin et al. 201 1), 
and more similar to the value found in the de nse gas surrou nding 
the UCHii region IRAS 20343-^4129 (.Fontani et al. 11201 2h . 



5.2. Comparison witli cliemical models 

In this section, we compare the observational results with the re- 
sults of different chemical models, in order to qualitatively dis- 
tinguish which effects may be dominant in the different phases 
that co-exist in this region. To fully model the chemistry of this 
complex is beyond the scope of this paper. Our goal is to explore 
the possible impact of different factors such as the grain-surface 
chemistry or time dependency, on the carbon chemistry in these 
environments. 

The PDR layers (hereafter, Lpdr and Lhd for the exposed 
PDR and the high density PDR, see Fig.O have different char- 
acteristics, in terms of physical conditions and time scales com- 
pared to the cold envelope. We assumed the density structure 
described in Sect.|2] (see also the top panel in Fig. (8]), in which 
the mainly atomic layer Lpdr (Ay < 1, = 2 x 10^ cm"^. 
Go = 5 X 10^) is followed by the high density layer Lhd 
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Table 4. Comparison of CCH and C-C3H2 abundances in different PDRs. 





Go 


[cm 3] 


X(CCH)(*) 


X(c-C3H2)(*) 


7V(c-C3H2)/A^(CCH) 


References 


MonR2 


MonR2 IF 






5.7 X 10-9 


3.5 X 10-11 


0.006 


a,b 


MonR2 MP2 






2.6 X 10-9 


7.6 X 10-11 


0.029 


a,b 


MonR2 FDR 


5x 10^ 


3x 10^ 


1.7 X 10-^ 


1.4 X 10-10 


0.008 


a 


MonR2 env 


1-100 


5x 10^ 


3.5 X 10-9 


1.4 X 10-11 


0.004 


a 


MonR2 red wing 






1.8 X 10-^ 


7.8 X 10-10 


0.04 


a 


Other FDRs 


Horsehead 


100 


1 X 10^ 


1.4 X 10-^ 


1.1 X 10-9 


0.08 


c 


Oph-W 


400 


2x 10^ 


6.0 X 10-9 


6.0 X 10-10 


0.1 


d 


1C63 


1500 


1 X 10^ 


1.1 X 10-^ 


9.0 X 10-10 


0.08 


d 


NGC7023 


2600 


10^ - 10^ 


6.0 X 10-9 


2.0 X 10-10 


0.03 


e,f,g 


Orion Bar 


2x 10^ 


1 X 10^ 


2.0 X 10-9 


3.8 X 10-11 


0.01-0.03 


f,h,i 



Refer ence s, (a): this work; (b) : iGinard et all (120121): (c): iFetv et all (l2005h (d) iTeyssier et all (I2004D : (e) iFuente et all (Il993h : (f): iFuente et all 
(l2003h: (g)lFilleri et al.l(i2012bh: fh):'Unger echts e75] (119971) : (i): IFuente et al]([l996) 

The uncertainty in the abundances is drive n by the uncer t ainty in the H2 gas column density, which is assumed to be a factor of 2. 
Calculated using the Ci^O observations of Gi nard et al.l (2012) as reference for the H2 gas column density: N(C^^O) = 1.3 x lOi^ cm-^. 
In MP2, we have to rely only on the 3mm and 1mm transitions of CCH. Thus, we may underestimate the CCH column density (and 
consequently, overestimate A/^(c-C3H2)/A/^(CCH)) by a factor of 2. 




N(c-C3H2)/N(CCH) 



Fig. 7. Color image shows the ratio of C-C3H2 to CCH column 
densities as derived from 1mm observations, assuming the rota- 
tional temperature of Fig.|6l Pixels without a detection of both 
CCH and C-C3H2 have been masked. The black contours repre- 
sent the integrated intensity of the 262.004 GHz CCH line (from 
15 to 25 Kkms"^ in steps of 2.5 Kkms"^), and white contours 
represent the Spitzer IRAC 8 fim emission. 



(1 < Ay < 10, nii2 = 3 X 10^ cm"^), where many molec- 
ular species start to be abundant. These innermost layers are 
surrounded by the low-density envelope Lenv (10 < Ay < 50, 
=5x10"^ cm"^). The assumption of constant density has to 
be taken with caution: local variations of gas density (clumps) 
and smoother gradients are expected to be found in the FDR lay- 
ers, and the density in the envelope most likely decreases with 
the distance from the IF. However, this zero-order approximation 
allows to reduce the free parameters of the modeling to a few, 
which ease the comparison of models and observational data. 



To mo del the LppR and Lhd layers, we used the Meudon 
FDR code dLe Petit et a l. 2006; Goicoechea & Le Bour loil2007 l : 
iGonzalez Garcia et al.l I2OO81: iLe Bourlot et al.i i2012) version 
1.4.3 using the density profile shown in Fig.[8l The Meudon 
Code computes the steady- state solution to the thermal balance 
and gas-phase chemical networks using accurate radiative trans- 
fer calculations and a plane-parallel geometry. The transition be- 
tween Lhd and Lpdr is calculated self-consistently by the code. 
This yields stationary solution for the gas and dust tempera- 
ture, as well as the abundance of each molecule in each of the 
layers by assuming equilibrium conditions. Time-dependent ef- 
fects and gas-grain interactions are not included in the model 
and therefore its use is justified for the region where chemical 
equilibrium has already been reached and where gas-grain inter- 
action is not expected to play a major role. However, the Meudon 
FDR code has very extensive gas-phase chemical networks and 
is best-suited for studying the gas-phase chemistry of Mon R2. 

To study the l ow-density, col d envelope we used the 
UCL.CHEM code dViti et al.1 l2004l) . This code includes gas- 
grain interaction: for the purpose of this paper it is used as a two- 
phase model which consists of the collapse of a prestellar core 
(Phase I), followed by the subsequent warming and evaporation 
of grain mantles due to the enhanced (with respect to the ISRF) 
radiation field (Phase II). In phase I, a diff'use cloud of density 
10^ cm"^ undergoes free-fall collapse until it has reached a fi- 
nal density (set by the user). This phase occurs at a temperature 
of 10 K. During the collapse, atoms and molecules collide with, 
and freeze onto, grain surfaces. The advantage of this approach 
is that the ice composition is not assumed, but it is derived by 
a time-dependent computation of the chemical evolution of the 
gas-dust interaction process which, in turns, depends on the den- 
sity of the gas. Hydrogenation occurs rapidly on grain surfaces, 
so that, for example, some percentage of carbon atoms accreting 
will rapidly become frozen out as methane, CH4. 

For all the models we use the same initial elemental abun- 
dances (see Table [5]) . The codes employ the reaction rate data 
from the UMIST astrochemical database, augmented with g rain- 
surface (hydrogenation) reactions as in IViti et al.1 (l2004l) and 
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IViti et alJ (l2Qllh for the UCL.CHEM models. In both Phases 
of this code, non-thermal desorption is also considered as in 
[Roberts et aLl(l2QQ7h . 



Table 5. Input parameters for the Meudon FDR code 







Value 




Radiation field intensity 


J X iU 


G- 


External Radiation field intensity 


100 


Av 


Total cloud depth 


50 mag 


Extinction 




Standard Galactic 






3.1 




Cosmic ray ionization rate 


5 X 10-1^ s-i 




dust minimum radius 


3 X 10-^ cm 


^max 


dust maximum radius 


3 X 10-^ cm 


a 


MRN dust size distribution index 


3.5 


He/H 


Helium abundance 


0.1 


0/H 


Oxygen abundance 


3.2 X 10-^ 


C/H 


Carbon abundance 


1.3 X 10-4 


N/H 


Nitrogen abundance 


7.5 X 10-5 


S/H 


Sulfur abundance 


1.8 X 10-5 




10- 



10-^" 
10^ 

o 

2 10° 
>^ 
'c75 

^ 10-' 

o 

O 













C,H 



1 


CH/CCH 

u — — 












C-C3H2/CCH 







f 



5.2.1 . Abundances in the PDRs 



0.0810 

Physical distance [pc] 



The column densities and abundances in the PDR obtained from 
the observations refer to the innermost 10 mag of the PDR. 
Therefore, we have to compare them with the mean values in 
the LpDR -\-Lhd layer. Another way to study the chemical differ- 
ence between these species and minimize the uncertainties due 
to the assumed physical structure is to consider the column den- 
sity ratios A^(CH)/A^(CCH) and A^(c-C3H2)/A^(CCH) instead of 
absolute abundances. 

In the exposed PDR layer (Lpdr), the gas temperature and 
UV field are so high that gas-phase chemistry is expected to 
be dominant. We can safely assume that all hydrocarbons have 
evaporated from grain surfaces an d that there are no i ces, since 
the dust temperature (100-150 K, iPilleri et al.ll2Q12ah is above 
the sublimation limit for all the molecules involved in reactions 
with hydrocarbons, especially methane (Tielens 2005) . 

The results of the Meudon code modeling are displayed in 
Fig.[8l The model shows that all the hydrocarbons are expected 
to have their abundance peak at Ay ~ 1. For CCH, the model 
predicts a peak abundance of ~ 10"^ which rapidly decreases to 
10"^^. Similarly, C-C3H2 has an abundance peak toward the PDR 
of ~ 10"^^ and a decrease to 10"^^. The abundance gradient is 
less abrupt for CH which passes from a peak of 10"^ in the PDR 
to 10"^^ in the more shielded layers of the molecular cloud. 

The column densities (and abundances) of CH and CCH are 
relatively well reproduced by the Meudon PDR code, consider- 
ing the integrated column density in the two layers Lpdr and 
Lhd (see Table O. Concerning C-C3H2, the model predicts a 
A^(c-C3H2)/A^(CCH) ratio of 0.0003 while the observed value is 
of the order of 0.01 in the PDR, which shows that the C-C3H2 is 
overabundant compared to the model predictions. This is consis- 
tent with several previous studies showing that in PDRs the gas- 
phase model predictions canno t account for the observed abun- 
dances of small h ydrocarbons (iFuente et al.ll 19931: iTeyssier et al.l 
120041: lPetvet"aDr20051 especially C-C3H2. 



Fig. 8. Model results using the gas-phase steady- state chemistry 
with the Meudon PDR code. The cloud is illuminated from the 
left side with a UV field of Go = 5 x 10^ 

A link with PAHs? 

A possible way to reduce the discrepancy could be the injection 
of fresh hydrocarbons by the photo-destruction of PAHs or larger 
grains (Fuente et al. 2003; Pety 2005). Indeed, laboratory exper- 
iments show that the photodissociation of small PAHs (A^c ^ 24) 
can lead to the production of small hy drocarbons C^H^ (n > m) 
that can be injected into the gas phase (lUseli Bacchitta & JoblinI 
2007). To understand if PAHs can be linked to a chemical difiTer- 
entiation of CCH and C-C3H2, it is useful to compare the ratio 
between the column density of C-C3H2 and CCH with the spatial 
distribution of PAHs. In Fig.[7] we display the comparison of the 
A^(c-C3H2)/A^(CCH) ratio with the IRAC 8 jim map. 

In low- to mild-UV irradiated PDR, the 8 jim map is dom- 
inated by the emission of the AIB at 1.1 [im, but in high-UV 
excited environments larger grains at thermal equilibrium can 
reach temperatures high enough to have continuum emission at 
these wavelength. To check whether this is the case in Mon R2, 
we have compared the IRAC 8 yum map with that obtained with 
the IRAC 3.5 yum filter, which is expected to be dominated ei- 
ther by the PAH 3.3 //m feature or by the stellar continuum. The 
two IRAC maps are very similar, except for very point-like fea- 
tures in correspondence with the infrared source IRS 2 and IRS 
3. Therefore, we can assume that most of the IRAC 8 jim emis- 
sion is indeed due to AIB emission. Finally, visual inspection 
of the IRS spectra in the central region of Mon R2 shows that 
the PAH features dominate the emission at 8 yum, except towards 
the bright infrared sources. The comparison of the ratio of the 
C-C3H2 and CCH column densities with the spatial distribution 
of CCH suggests that /) C-C3H2 is overabundant in the inner- 
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most layers of the cloud, and //) the C-C3H2 spatial distribution 
extends to the NW, following the northern IR ridge that is ob- 
served at 8 yum ( G inardet all2012h . 

The correlation between the PAH emission and the C-C3H2 
overabundance, although suggestive of a chemical link between 
these species, does not prove yet that these two tracers are chem- 
ically linked (i.e. by UV processing), since they are both good 
FDR tracers. For instance, it has been observed th at PAH and 
[C n1 emission are spatially correlated in PDRs (iJoblin et al.l 
l20ld) . since both tracers arise in the most exposed layers of the 
PDR. This could also mean that the chemical differentiation of 
CCH and C-C3H2 in PDRs can also be due to an increased abun- 
dance of C^ in these environments. Furthermore, C 11 and PAHs 
also trace the electron density, which is mainly due to the ioniza- 
tion of C and photo-electric effect on PAHs. Since C-C3H2 is sup- 
posed to form by electron recombination of the C3H3 , a higher 
production efficiency of photo-electrons ejected from PAHs may 
be a cause of C-C3H2 enhancement. However, the inclusion of 
these effects in PDR models is still very challenging and beyond 
the scope of this paper. 



one order of magnitude: for ^ > 10^ years, the agreement of 
A/^(c-C3H2)/A/^(CCH) with observations improves, but the abso- 
lute abundances and the CH/CCH ratio do not agree as well as 
at 1 X 10^ years. 

Repeating the calculations for higher freeze-out fractions 
(25% and 45%, see Table© one can find a better agreement 
for one or the other molecule, but it is difficult to find a solution 
that fits reasonably well all the abundances or column densities. 

In any case, the chemistry of Mon R2 seems to be dominated 
by grain-surface processes and time-dependent effects, since at 
the age of UCHII region (~0.1 Myr) the molecular chemistry 
is far from the steady state yet. In this scenario the current 
abundances of small hydrocarbons are dependent on the previ- 
ous history of the cloud. This is consistent with recent results 
of Mook erjea et al.l (12012*), who studied these effects in the en- 
velope associated with DR21. They investigated the impact of 
time-dependent effects for the carbon chemistry of this warm 
cloud, showing that PAH-related chemistry is not needed to ex- 
plain the abundances of neither CCH or C-C3H2 in this molecular 
envelope. 



5.2.2. Low density envelope 

Our observations indicate that the abundances of small hydrocar- 
bons in the envelope are about 5-10 times lower compared to the 
PDR. The column densities of CH and C-C3H2 relative to CCH 
are also a factor of 3 and 2 lower compared to the PDR. Even 
though, the CCH and C-C3H2 abundances are larger by several 
orders of magnitude than those predicted by steady-stage gas 
chemical models. In our sketch of MonR2, we assume the enve- 
lope to be 40 mag thick on each side of the spherical cloud, with 
a gas density of =5x10^ cm"^ and a gas kinetic temperature 
of 35 K. With these physical conditions, gas-grain chemistry and 
time-dependent effects may play a role. 

To study these processes we used the UCL_CHEM models 
running Phase 1 from diffuse core up to a density of = 
5x10^ cm"^, while the temperature and radiation field were kept 
at lOK and 1 Habing respectively. Once the final density was 
reached we increased the radiation field to Go = 100 and the ki- 
netic temperature to 35 K (Phase 2). The model assumes a num- 
ber of parameters, such as the freeze-out, photo- and thermal- 
desorption efficiencies. We use the standard valu es for the dif- 
ferent non thermal desorption efficiencies as in [Roberts et al.l 
(I2OO7I) . In particular, the models are strongly dependent on the 
freeze-out efficiency during phase 1 , which determines the ini- 
tial abundances of CH, CCH and C-C3H2 in Phase 2. We have 
tested the impact of the kinetic temperature in the envelope in 
the final hydrocarbon abundances, and found no significant dif- 
ference in the range [10-35] K. This is consistent with the fact 
that at rdust=35 K, m ost of the species remain in the grain man- 
tles (IViti et al.ll2004l) and the rates of main gas phase reactions 
involved in the formation and destruction of hydrocarbons do 
not present strong variations in this temperature range. 

In Fig. [9] we report an example of the time evolution during 
Phase 2 for an Ay = 41, assuming a low (10%) CO freeze- 
out. The model predicts that all the hydrocarbon abundances 
are constant up to 10"^ years. At this time, each of the abso- 
lute abundances for CH and CCH are reproduced within a fac- 
tor of 3 but their ratio is over-predicted by a factor of 7 (see 
Table [6]). Considering the high uncertainties in the absolute col- 
umn density of CH (driven by the unknown excitation temper- 
ature) it is difficult to further constrain the model based on this 
ratio. After 1 x 10"^ years, gas phase CH decreases by a factor 
of a few, whereas CCH and C-C3H2 decrease by the more than 



Av = 41 mag 



X(CCH)g,3 








X(CCH)g,,,3 X(CH)g,3 

X(C3H2)gas 




"X(C3H,V,,3X10^ 



X(CH)/X(CCH) 








X(C3H2)/X(CCH) 





Time [yrs] 



Fig. 9. Results of the UCL CHEM code to model the low-density 
envelope Lhd for a starting CO freeze-out of 10%. The results 
are given at Ay of 41 magnitudes and by setting the temperature 
fixed to 35 K. The bottom panel shows the ratio of gas-phase CH 
and C-C3H2 relative to CCH. The vertical dashed line marks 
the age of 10"^ years, which is used in Table [6] 



6. Conclusions 

In this paper, we have presented new observations of small hy- 
drocarbon molecules (CH, CCH and C-C3H2) in MonR2 ob- 
tained with the IRAM-30m telescope and the HIFI instrument 
onboard Herschel. In particular, the HIFI observations allowed 
to access the N=6-5 rotational transitions of CCH with a good 
spatial and spectral resolution. This allowed to constrain the 
excitation conditions and to distinguish the different phases in 
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Table 6. Comparison between the abundances derived from observations and model predictions. For the UCL_CHEM code, we 
show the results at ^ = 10"^ years. 







Ay 

[mag] 


X(CH) 


X(CCH) 


X(C-C3H2) 


A^(CH)/A^(CCH) 


A^(c-C3H2)/A^(CCH) 


vPDR 

^ohs 




1-10 


3.3:6-0 X 10-9 


1.7 X 10-^ 


1.4 X 10-10 


0.15 


0.008 


Meudon Lpor +Lhd 




1-10 


1.1 X 10"^ 


4.8 X 10"^ 


1.3 X 10"^^ 


0.22 


0.00025 


^env 
obs 

UCL_CHEM Lenv f = 
UCL.CHEM Lenv f = 

UCL_CHEML,„,f = 


lO(^) 
45(^) 


10-50 
41 
41 
41 


\.9^\l X 10-10 
5.1 X 10-10 
9.5 X 10-10 
7.5 X 10-10 


3.5 X 10-9 
1.3 X 10-9 
6.9 X 10-10 
3.5 X 10-10 


1.4 X 10-11 
4.8 X 10-13 

3.5 X 10-11 
1.4 X 10-10 


0.054 
0.39 
1.37 
2.1 


0.004 
0.0004 
0.05 
0.40 



Notes. For the observed abundances and the Meudon PDR results, we report mean values of the abundances, while for the UCL_CHEM results, 
we report the value at a given Ay, and assume that these values are relatively constant in the respective layers. 

Except for CH, where the uncertainty is driven by the excitation temperature, the abundance uncertainty is driven by the H2 gas reference, and 
is a factor of 2. 

The / value indicates the percentage of CO on the icy mantles at the end of Phase I. 



which this species is abundant. Profiting of a multi-transition and 
spatial study, we have derived the spatial distribution of these 
molecules, determining their abundances in the diff'erent layers 
of Mon R2. 

We compared the observational results with predictions of 
diff'erent PDR models, to investigate which kind of chemistry 
can dominate in the diff'erent environments that co-exist in 
MonR2: 

- In the high-UV illuminated PDR, grain mantles have already 
been evaporated and we can safely assume that chemical 
equilibrium has been reached. In this case, the gas-phase 
chemistry can explain quite well the abundances of CH and 
CCH but fails to reproduce C-C3H2. The chemical diff'eren- 
tiation of these two molecules happens in correspondence 
with PAH emission, suggesting once again a link between 
these species. 

- In the cold (~ 35 K) low-density molecular cloud, time 
dependent eff'ects seems to be very important: around 10"^ 
years, the abundances predicted by the models are relatively 
consistent with those derived from the observations. At this 
time all the absolute abundances are relatively well repro- 
duced, although CH is over-predicted relative to both CCH 
and C-C3H2. These models are however strongly dependent 
on the initial abundances due to the freeze-out in the cloud 
collapse phase. 

This study shows that several eff'ects must be taken into ac- 
count to understand the chemistry of small hydrocarbons such 
as CCH and C-C3H2: the relative contribution of UV radiation, 
grain-surface chemistry and time-dependent eff'ects may vary in 
diff'erent environments. Regarding PDRs, a more detailed anal- 
ysis of the small hydrocarbon chemistry in these environments 
requires: (/) a larger sample of PDRs spanning a large range of 
physical conditions to minimize effects due to a peculiar geom- 
etry of a particular source; (//) to revise the possible chemical 
routes to the formation of small hydrocarbons such as the photo- 
destruction of PAHs, and very small grains. Herschel data will 
allow to constrain the hydrocarbons chemical models, includ- 
ing related species that are not observable from the ground (e.g., 
C^, CH^). Finally (///) higher resolution observations with the 
Plateau de Bure Interferometer, NOEMA, ALMA and JCMT 
will allow to observe several key species with an unprecedented 
spatial resolution. This will enable to resolve the small-scale 



chemical frontiers in the PDR that seems to be key interpret the 
link between small hydrocarbons and PAHs. 
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Appendix A: Diffuse gas in the line of sight 

The low-lying transitions of all the hydrocarbons show a clear 
self-absorption at about the velocity of -llkms"^ toward the 
IF (see Fig. [3]). A similar, narrow self- absorption compone nt has 
also been reported for the [C 11] 158 yum line ( Pilleri et al.ll2Q12al; 
lOssenkopf et al.ll2Q13l) . This feature suggests the presence of a 
relative diffuse phase in the line of sight. The velocity at which 
the absorption features appears is the same as that of the emis- 
sion of the quiescent molecular cloud. Thus, it is likely that this 
absorbing layer is a low-density gas associated with the most 
external layer of the molecular cloud. 

To determine the column density of the absorbing gas, we 
need to estimate a reference for the spectral profile of the emis- 
sion line. We assume that the widths of the lines belonging 
to a species are similar, and use a transition that is less af- 
fected by self-absorption to estimate the FWHM of a gaus- 
sian profile. We obtain a FWHM of 5.5 km s"^ for CCH (us- 
ing the J = 87.407 GHz line) and 3kms"^ for C-C3H2 (using 
the 150.820 GHz line) and 6kms-i for CH. Using these line 
widths, we can adjust the blue-shifted wing of the self-absorbed 
line with a gaussian and use it as reference to calculate the 
integrated opacity of the self- absorption. By assuming an ex- 
citation temperature Tex = 2.73 K, we derive column densi- 
ties associated to the foreground of N(CU) = 4.8 x 10^^ cm"^ 
A^(CCH) = 5.0 X 10^^ cm-2 and A^(c-C3H2) = 1.4 x 10^^ ^^^-2 

From CH, we can derive the H2 column density associ- 
ated with the absorption, hence an estimate of the extinction. 
The standa rd abundance of CH relative to H2 in diff'use gas is 
3.5 X 10"^ (ISheff^er et al.l l2008'). which leads to a molecular hy- 
drogen column density of A^(H2) ~ 1.3 x 10^^ cm"^ or about 
1.3 mag of extinction. This yields CCH and C-C3H2 abundances 
(relative to H2) of 3.9 x 10"^ and 1.6 x 10"^, respectively. These 
values a re consistent with the typical values o bserved in dif- 
fuse gas dLiszt & Lucasll2002l:lGerin et al.ll201lh . The results are 
similar when considering an offset of [20'',20''], where the self- 
absorption is still detected with a sufliicient S/N to infer relatively 
accurate column densities. 

These estimates suff'er from the uncertainty due to the as- 
sumed line profile and on the observed continuum. Nevertheless, 
the detection of this foreground, lAy thick layer associated with 
the external layer of the envelope is interesting and moves to- 
ward a deeper understanding of this source. 



15 



R Pilleri, et al.: Spatial distribution of CCH and C-C3H2 in the Ultra Compact HII region Monoceros R2 

Appendix B: Spectroscopic parameters, maps and 
rotational diagrams 



Table B.l. Spectroscopic parameters of the observed transitions 



Transition 




Frequency 


Em 


Ki 


gu 






[GHz] 


[K] 


[s-M 




CH 2Ui/2, N =l,J = 


3/2- 


1/2 








P = 2- - 1+ 




536.76115 


0.00072 


0.0638 


5 


P= 1- - 1+ 




536.78195 


0.00072 


0.0213 


3 


p = r - 0+ 




536.79568 


0.0 


0.0425 


3 


CH 2ni/2, N = 2,J = 


5/2- 


3/2 








P = 3+ - 2- 




1656.961 


25.727 


0.03719 


7 


P = 2+ - 2- 




1656.970 


25.727 


0.00372 


5 


P = 2+ - 1- 




1656.972 


25.727 


0.03348 


5 


CCH N = I - 












J=3/2-l/2, F= 1-1 




87.28410 


0.00216 


2.5920 X 10"^ 


3 


J=3/2-l/2, F= 2-1 




87.31690 


0.00216 


1.5274 X 10-6 


5 






/ .jz^ojy 


00000 


1 9^^84. V 1 0~6 


3 


J= 1/2- 1/2, F= 1-1 




87.40199 


0.00216 


1.2715 X 10-6 


3 


J= 1/2- 1/2, F= 0-1 




87.40716 


0.00216 


1.5324 X 10-6 


1 


CCH = 3 - 2 












J=7/2-5/2, F=4-3 




262.00426 


12.57513 


5.30542 X 10-5 


9 


J=7/2-5/2, F=3-2 




262.00648 


12.57383 


5.10470 X 10-5 


7 


J=5/2-3/2, F=3-2 




262.06499 


12.58203 


4.87506 X 10-5 


7 


J=5/2-3/2, F-2-1 




262.06747 


12.58261 


4.45975 X 10-5 


5 


CCHA^ = 6-5 












1=13/2-11/2, F=7-6 




523.97157 


62.87103 


4.5697 X 10-4 


15 


1=13/2-11/2, F=6-5 




523.97217 


62.86987 


4.5179 X 10-4 


13 


J= 11/2-9/2, F=6-5 




524.03390 


62.88685 


4.4939 X 10-4 


13 


J = ii/Z-9/Z, r=j-4 




CO /I /ICQ 


OZ.OO /J 1 


A /to AC -.y 1 rv— 4 
4.4zUj X iU 


1 1 


C-C3H2 Jka,Kc = 2i,2 


-lo,i 


85.33889 


4.1 


2.55335 X 10-5 


15 


C-C3H2 Jka,Kc = 4o,4 


-3i,3 


150.82067 


19.3 


1.79933 X 10-4 


9 


C-C3H2 Jka,Kq = 4i^4 


- 3 0,3 


150.85191 


17.0 


1.80086 X 10-4 


27 


C-C3H2 / = 6 - 5 












Ka,Kc = 1,6-0,5 




217.82215 


36.3 


5.93376 X 10-4 


39 


Ka,Kc =0,6-1,5 




217.82215 


38.6 


5.93403 X 10-4 


13 


C-C3H2 Jka,Kc - 5l,4 


-42,3 


217.94005 


33.1 


4.42680 X 10-4 


33 


C-C3H2 Jka,Kc = 43,2 


-32,1 


227.16913 


26.7 


3.42484 X 10-4 


27 


C-C3H2 / = 7 - 6 












Ka,Kc = 0,7-1,6 




251.31434 


48.3 


9.34953 X 10-4 


45 


Ka,Kc = 1,7-0,6 




251.31434 


50.7 


9.34781 X 10-4 


15 


C-C3H2 / = 6 - 5 












Ka,Kc = 2,5 -1,4 




251.52730 


45.1 


7.42318 X 10-4 


39 



References: CH: (McCarthy et al. 2006) & CDMS; CCH: 
(iPadovani et al.ll2009l) & CDMS: C-C3H2: JPL 
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Fig. B.l. Integrated intensity of all the observed C-C3H2 and CCH transitions. The upper state energy level and the principal quantum 
numbers are shown for each transition. The smaller-scale maps have been convolved to the spatial resolution of 15"to improve the 
S/N ratio. 
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